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ABSTRACT 


This  thesis  describes  the  development  of  a  system  for 
measuring  frequency-domain  scattered  fields  in  the  Transient 
Electromagnetic  Scattering  Range  at  the  Naval  Postgraduate 
School.  The  new  system  employs  a  stepped- frequency  CW 
waveform  and  utilizes  an  HP-8510B  network  analyzer  as  an  RF 
front-end  and  a  coherent  receiver.  A  pair  of  AEL  H1498 
antennas  was  installed  to  cover  a  frequency  range  of  2  GHz  to 
18  GHz.  An  HP-82300C  BASIC  Language  Processor  was  installed 
on  a  COMPAQ  Deskpro-386  PC,  and  an  HP-BASIC  program  was 
developed  for  remote  control  of  the  HP-8510B  with  data 
acquisition  over  the  HPIB  bus.  A  post-processing  algorithm 
was  created  using  MatLab  for  background  subtraction, 
calibration,  and  deconvolution.  A  set  of  RCS  measurements 
was  made  using  various  size  spheres,  and  the  post-processing 
outputs  were  compared  to  computed  values.  Good  agreement 
between  these  measurements  and  computed  data  indicates 
excellent  accuracy  of  the  measurement  system  and  valid 
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INTRODOCTIOM 


I  . 


A .  BACKGROUND 

A  radar  detects  targets  at  far  distances  by  transmitting 
radio  waves  and  listening  for  echoes  resulting  from  fields 
scattered  by  the  target.  The  scattered  fields  vary 
significantly  with  the  radar  frequency  and  polarization  and 
with  the  target's  shape,  materials,  and  orientation. 
Knowledge  of  the  scattering  properties  of  potential  targets 
is  essential  in  designing  and  assessing  performance  of  radar 
systems.  On  the  other  hand,  scattering  aspects  of  new 
military  vehicle  designs  must  be  considered  to  enhance  their 
survivability  against  enemy  radars  and  weapons.  Because 
radars  play  increasingly  important  roles  in  modern  warfare, 
research  to  develop  efficient  and  accurate  techniques  of 
measuring  scattering  properties  of  radar  targets  has  become  a 
high-priority  task. 

Given  a  knowledge  of  the  incident  wavefield,  the 
electromagnetic  properties  of  a  radar  target  are  completely 
described  by  the  scattered  wavefield.  The  strength  of  the 
scattered  field  is  indicated  by  the  target’s  radar  cross 
section  (RCS),  a  parameter  widely  used  in  the  radar 
community.  Although  the  RCS  is  a  normalized  scalar  measure 
of  the  complex  scattered  wavefield,  the  term  is  often  used  to 
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describe  the  scattering  characteristics  of  a  target. 
Scattering  responses  of  the  target  can  be  measured  either  in 
the  time  or  frequency  domain  and  the  two  responses  are 
related  by  a  Fourier  transform. 

The  Transient  Electromagnetic  Scattering  Laboratory 
(TESL)  at  the  Naval  Postgraduate  School  (NPS)  currently 
utilizes  a  time  domain  technique  to  measure  scattered  fields. 
The  system  transmits  short  pulses  and  measures  the  target 
scattering  response  directly  using  a  sampling  oscilloscope. 
Time-domain  techniques  are  considered  faster,  less  expensive, 
and  better  suited  to  direct  observations  of  transient 
scattering  responses  than  frequency-domain  techniques. 
However,  frequency-domain  techniques  are  considered  more 
stable,  less  susceptible  to  noise,  and  capable  of  spanning  a 
wider  measurement  frequency  range. 

B .  OBJECTIVE 

The  objective  in  this  thesis  research  is  to  develop  a 
capability  for  measuring  frequency-domain  scattered  fields 
over  the  2  to  18  GHz  range.  The  effort  consists  of 
installing  and  integrating  an  HP-8510B  network  analyzer 
system,  mounting  antennas  on  a  wall  panel  in  an  anechoic 
chamber,  installing  an  HPIB  controller  in  a  personal 
computer,  writing  an  HP  BASIC  program  for  remote  system 
control  and  data  acquisition,  creating  a  MatLab  signal 
processing  algorithm  for  background  subtraction,  calibration 
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and  transformation  to  the  time  domain,  and  conducting  RCS 
measurements  on  simple  targets  for  validation.  The  new 
frequency-domain  system  will  complement  the  existing  time- 
domain  RCS  measurement  system  and  should  yield  similar 
results  when  frequency -domain  data  are  transformed  into  the 
time  domain.  Measurement  validations  for  complex  targets 
will  be  possible  through  comparing  time  and  frequency-domain 
results  acquired  independently. 

The  thesis  is  divided  into  six  chapters-.  Chapter  ll 
illustrates  the  system  and  measurement  setup.  Detailed 
descriptions  are  given  for  the  anechoic  chamber,  the  HP-8510B 
network  analyzer,  and  the  cintennas.  Chapter  III  explains  the 
system  control  and  data  acquisition  including  controller 
hardware,  data  acquisition  software,  and  measurement 
procedures.  Chapter  IV  describes  signal  processing  software 
for  background  subtraction  and  calibration.  Fast  Fourier 
Transform  (FFT)  techniques  to  recover  time  transient 
responses  from  frequency  domain  data  are  discussed,  in 
addition  to  important  processing  aspects  of  windows, 
conjugation,  and  zero  padding.  Chapter  V  addresses  system 
validation  by  presenting  results  of  measurements  on  different 
size  spheres.  Finally,  conclusions  and  future  considerations 
are  presented  in  Chapter  VI. 


II.  SYSTEM  DESCRIPTION  AMD  MEASUREMENT  SETUP 


A .  OVERVIEW 

The  frequency-domain  scattered  field  measurement  system, 
depicted  in  Figure  1,  is  implemented  on  an  HP-8510B  network 
analyzer  [Ref.  1]  using  its  stepped  frequency  mode.  The 
system  consists  of  an  HP-8340B  synthesized  sweeper  as  an  RF 
source,  an  HP-8511A  frequency  converter  as  an  RF  front  end, 
an  HP-8510B  network  analyzer  as  an  IF  receiver,  AEL  H1498 
multi-octave  horns  serving  as  transmitting  and  receiving 
antennas,  and  an  IBM  PC-compatible  computer  with  HP-82300C 
BASIC  Language  Processor  as  a  data  acquisition  workstation 
and  remote  system  controller. 

System  components  are  selected  for  scattered  field 
measurements  of  2  to  18  GHz.  The  system  operates  by  spanning 
the  frequency  band  with  801  steps.  The  maximum  number  of 
frequency  steps  allowed  by  the  analyzer  is  used  to  maximize 
the  unambiguous  time  extent. 

For  each  frequency  step,  the  system  is  phase-locked  and 
continuous -wave  (CW)  signals  are  transmitted  through  the 
transmitting  antenna  in  an  anechoic  chamber.  The  RF  energy 
reflected  by  the  target  is  collected  by  the  receiving 
antenna,  and  initially  down-converted  to  20  MHz  using  a 
harmonic  mixer  in  the  frequency  converter,  followed  by  a 
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Computer 


Figure  1.  Schematic  Drawing  of  the  Frequency  Domain  Scattered  Field 
Measurement  System  Implemented  on  the  HP-8510  Network 


second  mixer  which  produces  an  IF  frequency  of  100  KHz  for 
subsequent  analcg-to-digital  (A/D)  conversion.  This  IF 
signal  is  processed  through  the  "bl"  channel  in  the  frequency 
converter.  For  coherent  detection  with  phase  information,  a 
constant  RF  signal  is  coupled  using  a  22  dB  directional 
coupler  at  the  RF  source  and  down-converted  through  the  “al“ 
channel.  Amplitude  and  phase  information  for  the  specific 
frequency  are  established  by  comparing  the  "bl"  channel 
(signal  channel)  to  the  “al“  channel  (reference  channel) . 
The  HP-8510B  periodically  sends  a  command  to  the  synthesizer 
to  change  to  the  next  frequency  step  until  all  801 
frequencies  are  measured. 

The  magnitude  and  phase  data  for  801  frequencies  are 
transferred  from  the  HP-8510B  to  the  computer  over  the  HPIB 
bus,  stored  on  the  internal  hard  dislc,  and  subsequently 
processed  with  a  digital  signal  processing  algorithm  written 
in  MatLab. 

B .  SYSTEM  DESCRIPTION 

1 .  Anechoic  Chamber 

The  anechoic  chamber,  shown  in  Figure  2,  is  located 
in  room  535  of  Spanagal  Hall,  measures  3.1  x  3.1  x  6.2  meter 
in  height,  width,  and  length  and  is  shielded  using  aluminum 
sheet  metal.  Antennas  mounted  on  the  source  wall  illximinate 
targets  supported  on  a  low-density  styrofoam  column 
positioned  2.3  meters  away.  The  source  wall  is  covered  with 


Metal  Enclosure 


Figure  2.  General  Layout  of  the  Anechoic  Chamber 

21  cm  pyramidal  radar  absorbing  material  (RAM)  and  the  back 
wall  is  covered  with  46  cm  pyramidal  RAM.  Sidewalls,  floor, 
and  ceiling  are  covered  with  longitudinal  wedge  RAM  to  direct 
energy  towards  the  absorbing  back  wall.  Figure  3  shows  the 
chamber  interior  looking  from  the  back  wall.  Further  details 
of  the  chamber  can  be  found  in  References  2  and  3. 

2 .  HP-8510B  Network  Analyzer 

The  HP-8510B  network  analyzer,  shown  in  Figure  4,  is 
a  high-performance  measurement  instrument  designed  to  make 
microwave  measurements  of  many  kinds.  In  the  course  of  this 
thesis  effort,  the  HP-8510B  network  analyzer  was  configured 
as  a  phase- locked  scattered  field  measurement  system  that 
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Figure  3.  Interior  of  the  Anechoic  Chamber 


Figure  4.  HP-8510  System  Mounted  on  a  Rack 
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offers  very  broad  frequency  coverage,  wide  <^TMuiiic  range,  and 
high  measurement  sensitivity.  The  HP-8510B  network  analyzer 
system  consists  of  a  microwave  source  (HP-83 40B) ,  RF  front 
end  (HP-8511A) ,  and  the  HP  8510B  as  IF  receiver  and  system 
controller.  The  HP-8510B  controls  the  RF  source  and  test  set 
over  a  dedicated  HP-8510B  system  bus.  The  HP-8511A  RF  front 
end  provides  a  broad  bandwidth  of  45  MHz  to  26.5  GHz  and  wide 
dynamic  range  of  75  to  105  dB,  depending  on  averaging  factor. 
The  HP-8511A  uses  a  harmonic  sampling  technique  to  convert 
the  RF  signals  at  its  two  inputs  to  20  MHz  IF  signals  that 
are  subsequently  processed  by  the  HP-8510B. 

In  manual  operation,  the  front  panel  keys  of  the  HP- 
8510B  are  used  to  define  measurement  parameters  and  control 
the  measurement  process.  In  remote  operation,  the  personal 
computer  controls  the  HP-8510B  by  sending  commands  over  the 
HPIB  bus  to  execute  the  measurement  process. 

To  measure  amplitude  and  phase,  the  HP-8510B  must 
have  a  reference  signal  that  remains  constant  during  the 
measurement.  The  phase-locked  reference  signal  is  obtained 
using  an  HP-11691D  20  dB,  2-18  GHz  directional  coupler 
inserted  between  the  HP8340B  RF  source  and  the  transmitting 
antenna.  The  measured  amplitude  and  phase  data  are 
transferred  to  the  personal  computer  over  an  HPIB  bus. 

For  wideband  scattered  field  measurements,  the  HP- 
8510B  network  analyzer  can  operate  either  in  RAMP  mode  or  in 
STEP  mode.  In  RAMP  mode,  the  source  is  swept  in  a  continuous 
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analog  sweep  from  the  lower  to  upper  frequency  and  the  data 
is  sampled  without  stopping  the  sweep.  In  STEP  mode,  the 
source  is  tuned  and  phase-locked  to  each  of  the  801  frequency 
sample  points.  The  STEP  mode  is  more  stable  in  terms  of 
frequency  accuracy  auid  source  power  level  than  the  RAMP  mode. 
The  STEP  mode  will  be  used  in  all  results  of  this  thesis. 

3 .  Antennas 

AEL  H1498  multi-octave  horns  [Ref.  4]  are  used  for 
transmit  and  receive  antennas  to  meet  wideband  requirements. 
This  coaxially-fed,  do\ible-ridged  horn  antenna  covers  the  2- 
18  GHz  frequency  band.  Average  electrical  characteristics 
are  listed  in  Table  1  and  performance  over  the  2-18  GHz 
frequency  range  is  shown  in  Figxire  5. 


TABLE  1.  ELECTRICAL  CHARACTERISTICS  OF  AEL  H1498  ANTENNA 


Frequency 

(GHz) 

VSWR 

Avg.  Gain 
(dBi) 

Avg.  3dB  Bearawidth 

E  Plane  H  Plane 

WEiSSSSmM 

2-18 

HHHIi 

11 

o 

O 

— 

50 

The  antennas  are  mounted  on  a  24 -inch  square  wood 
panel  covered  with  aluminum  sheet  which  is  tightly  fitted 
into  a  square  opening  in  the  source  wall.  The  panel  is 
covered  with  the  saune  radar  aibsorbing  material  as  the  source 
wall  to  eliminate  multi -path  effects  between  the  antennas  and 
the  panel.  The  antennas  are  displaced  vertically  with  the  E- 
plane  in  the  horizontal  direction  to  minimize  coupling.  They 
are  shimmed  to  point  the  beam  at  the  target  to  minimize 
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B’igure  h.  Performance  Cvjrves  fo 


boresi^ht  errors .  Figure  6  shows  the  panel  covered  with  RAM 
material . 
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Figure  6.  Masking  of  Panel  with  Radar  Absorbing  Material 
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III. 


REMOTE  CONTROL  AMD  DATA  ACQUISITION 


A .  REMOTE  CONTROLLER  HARDWARE 

The  HP-8510B  is  used  to  define  measurement  parameters  and 
control  the  measurement  process;  it  can  be  operated  manually 
with  front  panel  keys  or  remotely  through  an  HPIB  bus.  In 
remote  operation,  measurement  parameters  such  as  frequency 
range,  number  of  frequency  points,  and  averaging  factor,  can 
be  set  from  a  computer  by  sending  proper  commands  on  the  HPIB 
bus.  A  COMPAQ  Deskpro-386  computer  using  an  HP-82300C  BASIC 
Language  Processor  [Ref.  5]  with  a  built-in  HPIB  interface 
was  chosen  as  a  remote  system  controller. 

The  HP-823 OOC  BASIC  Language  Processor  consists  of  an 
MC68000  processor  card,  HP-BASIC  interpreter,  emulator 
software,  4  megabytes  of  RAM  memory,  and  a  built-in  HPIB 
interface.  The  processor  card  fits  into  a  standard  16  bit 
ISA  slot  of  the  computer.  The  HP-BASIC  interpreter  and 
emulator  software  provide  the  capabilities  of  an  HP-300  BASIC 
workstation  while  retaining  all  the  capabilities  of  MS-DOS. 
Bi-directional  communication  between  HP-BASIC  and  DOS 
applications  allows  data  acquired  by  BASIC  to  be  stored  in 
DOS  ASCII  format  on  the  PC  hcurd  disk. 


15 


B .  CONTROL  AND  DATA  ACQUISITION  SOFTWARE 

Control  and  data  acquisition  software,  named  “GET_DATA“, 
is  written  in  HP-BASIC  and  runs  on  the  HP-82300C  BASIC 
Language  Processor  within  the  PC.  The  software  controls  the 
HP-8510B,  acquires  measured  data  through  the  HPIB  bus,  and 
writes  this  to  the  PC  hard  disk. 

After  the  PC  is  booted  up,  the  user  types  "cd  blp“  to 
change  to  the  blp  directory.  The  name  “blp*  stands  for  basic 
language  processor,  and  is  a  default  name  given  when  the 
BASIC  Language  Processor  is  installed.  The  user  then  types 
"basic"  to  cause  the  HP-BASIC  interpreter  to  run  on  the  PC, 
From  this  point,  the  PC  r\ins  as  an  HP-300  BASIC  workstation, 
and  the  user  issues  HP-BASIC  commands  until  quitting  HP-BASIC 
by  typing  " [control] -FI 0"  and  returning  to  the  DOS 
environment . 

When  HP-BASIC  sTuns,  the  user  types  "LOAD  'GET_data"". 
This  BASIC  command  loads  the  program  from  the  hard  disk  to 
RAM.  When  the  user  types  "RUN",  the  program  executes  and 
takes  over  control  of  the  HP-8510B.  The  preceding  executions 
are  sxammarized  as  follows: 

c:\>  cd  blp 

c:\BLP>  basic 

LOAD  "GET_DATA" 

RUN 
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When  the  “GET_DATA“  program  runs,  the  computer  CRT  displays 
the  main-menu  function  keys  at  the  bottom  and  the  current  HP- 
8510B  status  parameters  such  as  start  frequency,  stop 
frequency,  nxamber  of  frequency  points,  RF  source  power  level, 
sweep  mode,  active  calibration  set,  and  averaging  factor  at 
the  top.  The  user  can  recall  a  default  setup  and  start 
measurements  by  pressing  appropriate  function  keys. 


The  main-menu  function  keys  are  shown  in  Figure  7 ,  and 
their  fxinction  is  described  below. 


Figure  7.  Main-menu  Function  Keys 

Continu.! 

Sweep  -  Pressing  this  function  key  causes  the  HP-8510B 
to  sweep  continuously  through  the  801  frequency  points,  to 
transmit  RF  energy  and  update  data. 


-  Pressing  this  function  key  causes  the  HP-8510B  to 
start  sweeping  from  the  first  frequency  point  to  the  last 
frequency  point,  when  the  single  sweep  is  completed,  the  RF 
transmission  is  shut  off  and  the  data  are  stored  in  the  HP- 
8510B  until  the  "Continu.  Sweep"  key  or  "Single  Sweep"  key  is 
pressed  for  new  data. 

-  Pressing  this  function  key  halts  the  frequency 
sweeping  and  shuts  off  RF  transmission. 
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-  Pressing  this  function  key  alters  the  averaging 
used  to  increase  signal-to-noise  ratio.  In  STEP  mode, 
averaging  is  done  at  each  frequency  as  the  HP-8510B  steps 
through  the  801  frequency  points.  The  averaging  factor  can 
range  from  1  to  4096.  The  measurement  time  increases 
proportionally  as  the  averaging  factor  increases.  The  data 
in  the  results  section  were  measured  with  the  averaging 
factor  of  128. 

-  Pressing  this  function  key  causes  the  writing  of 
two  data  files.  The  first  file  begins  with  a  letter  "D"  and 
is  in  a  DOS  ASCII  format,  the  second  file  begins  with  a 
letter  “F"  and  is  in  an  HP-BASIC  binary  format.  The  contents 
of  these  two  files  are  exactly  the  same.  Each  line  of  801 
rows  lists  a  frequency  index,  magnitude  in  dB,  and  phase  in 
degrees.  After  two  plots  are  drawn  on  computer  CRT,  the 
program  is  paused  until  a  "continue"  function  key  is  pressed 
to  resume. 

-  Pressing  this  fxinction  key  allows  a  user  to 
retrieve  old  data  and  display  them  on  a  computer  CRT.  The 
program  reads  the  HP-BASIC  binary  format  file. 

-  Pressing  this  fxinction  k^  allows  the  user  to 
change  measurement  parameters  on  the  HP-8510B  remotely. 


-  Pressing  this  function  key  sets  measurement 
parameters  to  default  values.  The  default  values  are  2-18 
GHz,  801  frequency  points,  STEP  mode,  continuous  sweep,  a 
source  power  of  10  dBm,  and  an  averaging  factor  of  16. 


when  the  “Change  Setup"  function  key  is  pressed,  a  setup 
sub-menu  will  be  displayed.  The  function  keys  for  the  setup 
sub-menu  are  shown  in  Figure  8. 


Save 

State 


Recall  I  I Change  ||No.  of 
State  I IFreq  |  {Point 


Figure  8.  Setup  S\ib-menu  Function  Keys 


Save 

State 


-  Pressing  this  function  key  saves  a  complete 


HP-8510B  state  into  a  storage  register  in  the  HP-8510B.  All 
measurement  parameters  are  included  in  the  saved  state. 


Recall 

State 


-  Pressing  this  function  key  recalls  a  saved  state. 


Change 

Freq  -  Pressing  this  function  key  allows  the  user  to 

change  frequency  range. 


No.  oj 
Point 


-  Pressing  this  function  key  allows  the  user  to  set 


the  number  of  frequency  points  to  51,  101,  201,  401  or  801. 


Power 

Level 


-  Pressing  this  function  key  allows  the  user  to  set 


the  power  level  of  the  RF  source. 


[Sweep 

[Mode 


-  Pressing  this  function  key  allows  the  user  to 
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select  either  STEP  mode  or  RAMP  mode. 

-  Pressing  this  fxmction  key  activates  the  HP-8510B 
internal  error  correction  feature. 

-  Pressing  this  function  key  causes  the  main  menu 
to  be  displayed. 

The  program  can  be  stopped  by  pressing  ■ [shift] -pause" . 
The  source  code  may  be  edited  by  typing  "EDIT".  The  program 
must  be  stopped  before  "EDIT"  is  executed.  Further  details 
of  the  HP-BASIC  programming  can  be  found  in  Reference  1,  and 
the  "GET_DATA"  source  code  can  be  fovind  in  Appendix  A, 

C .  MEASUREMENT  PROCEDURES 

The  first  step  is  to  power  up  the  HP-8510B  system  and 
boot  up  the  coit5)uter  and  allow  a  30 -minute  warm-up  for  the 
synthesizer.  Next,  load  the  "GET_DATA"  program.  When  the 
program  rxins,  press  the  "Default  Setup"  function  key  to  set 
the  HP-8510B  to  a  default  setup.  Use  the  "Change  Setup"  key 
to  change  the  setup.  For  the  data  in  this  thesis,  the 
default  setup  is  first  recalled,  and  then  the  averaging 
factor  is  changed  to  128. 

when  the  setup  is  completed,  press  "Hold  Sweep"  key  to 
halt  RF  radiation.  This  step  may  be  skipped  if  the  user 
believes  the  RF  radiation  level  is  not  a  health  hazard. 
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All  measurements  require  three  steps:  calibration, 
background,  and  target  scattering.  In  the  first  step,  place 
a  metallic  sphere  calibration  target  on  the  target  column. 
Press  “Single  Sweep*  key.  Sweeping  through  all  801  points 
takes  40  seconds  to  200  seconds  depending  on  the  averaging 
factor.  When  sweeping  is  done,  press  “Record  Data*  key. 
When  two  plots  are  displayed,  press  “Continue*  key  and  type 
in  a  number  such  as  123  for  filenames  of  "D123.DAT*  and 
“F123.DAT*.  The  input  statement  is  programmed  to  accept  only 
numeric  inputs  and  not  alphanumerics .  If  the  user  desires 
other  names,  the  DOS  “rename*  command  can  be  used  in  the  DOS 
environment  after  the  data  acquisition  program  is  terminated. 
Remove  the  calibration  target,  and  measure  the  background  by 
pressing  the  “Single  Sweep*  key.  Press  the  “Record  Data*, 
and  assign  a  new  number  to  the  background  data.  Then,  place 
the  target  on  the  column  and  repeat  the  single  sweeping  and 
recording  for  the  target  data. 

Three  data  files  have  been  saved  on  the  PC  hard  disk; 
calibration,  background,  and  target.  These  files  will  be 
subsequently  processed  by  the  MatLab  program.  If  more  target 
measurements  are  required,  the  preceding  three  steps  can  be 
repeated  for  a  new  target,  using  new  names  for  the  stored 
files. 
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IV.  POST-PROCBSSINQ 

JL .  ALGORITHM  OVBRVIBW 

The  RCS  of  a  target  is  a  function  of  frequency, 
polarization,  and  target  orientation  and  should  be 
independent  of  measurement  instrumentation  and  anechoic 
chamber  conditions.  However,  the  raw  data  obtained  from  RCS 
measurements  contains  unwanted  effects  such  as  noise,  antenna 
coupling,  chamber  clutter,  antenna  gain,  and  distance  to 
tcirget . 

The  unwanted  effects  fall  into  two  groups;  the  first  is 
caused  by  extraneous  signals  such  as  noise,  antenna  coupling 
and  cheunber  clutter,  and  the  second  is  caused  by  the 
dependence  of  the  measurement  system  on  transmitting  power 
level,  cable  loss,  antenna  gain,  distance  to  target,  and 
receiver  mixer  performance.  In  an  ideal  case,  the  receiver 
output  in  the  absence  of  the  target  should  be  zero.  Due  to 
residual  scattering  of  chamber  RAM  and  significant  coupling 
between  the  transmitting  and  receiving  antennas,  however,  an 
undesirable  background  signal  exists  in  all  measured  data. 
This  undesirable  signal  can  be  removed  by  measuring  the 
background  data  and  performing  a  complex  subtraction  from  the 
target  data.  The  dependence  on  the  measurement  system  can  be 
removed  by  measuring  a  calibration  sphere  of  known  scattering 


22 


characteristics  and  the  target  of  interest  cmd  forming  the 
ratio  of  the  measured  responses  multiplied  by  the  computed 
field  response  of  the  calibration  sphere. 

Once  the  frequency  domain  scattered  fields  of  the  target 
are  established  by  the  background  subtraction  and  calibration 
process,  transient  time  responses  are  obtained  through  an 
inverse  Fourier  transform.  For  a  time  transient  response 
using  efficient  inverse  FFT,  con^lex  conjugates  and  zeros  are 
appended  to  the  frequency-domain  data,  as  will  be  shown. 
Finally,  time-gating  the  transient  response  around  the  target 
and  Fourier  transformation  back  to  the  frequency  domain 
smoothes  out  the  frequency  data  and  provides  refined  RCS 
levels . 

B .  BACKGROUND  SUBTRACTION  AND  CALIBRATION 

Operation  of  the  measurement  system  is  represented  by 
the  signal  flow  diagram  in  Figure  9.  The  RF  excitation  is 


Figure  9.  System  Signal  Flow  Diagram 
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denoted  ty  X,  euid  the  signal  out  of  the  network  analyzer  1:^ 
Y.  The  symbol,  Hs,  represents  the  transfer  function  of 
combined  measurement  system  dependence  which  includes  RF 
cable  loss,  both  antenna  characteristics,  two-way  free-space 
loss  from  antennas  to  target,  and  RF  mixer  conversion  loss. 
The  symbol,  Ht,  is  the  frequency  response  of  the  target,  and 
He  is  that  of  the  calibration  sphere.  Both  Ht  and  He  are  the 
ratios  of  the  scattered  field  at  the  location  of  the 
receiving  antenna  to  the  incident  field  at  the  location  of 
the  scatterer,  multiplied  by  the  one-way  free-space  spread 
factor,  as  defined  in  Equations  (1)  and  (2) . 


E®  (at  the  receiving  antenna) 

Ht  =  R  (4ji)i/2 - - - —  (1) 

E^  (at  the  target) 

E®  (at  the  receiving  antenna) 

He  =  R  (4jC)i/2  - - - - -  (2) 


(at  the  calibration  sphere) 

where  Ht  =  frequency  response  of  target 

He  =  frequency  response  of  calibration  sphere 

E®  =  scattered  electric  field 

E^  =  incident  electric  field 

R  =  distance  from  the  receiving  antenna  to 
the  target 

Background  errors  such  as  antenna  coupling  and  chamber 
clutter  are  grouped  together  and  denoted  by  Hb-  The  chamber 
thermal  noise  and  receiver  mixer  noise  are  denoted  by  N. 
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Equations  (3) -(5)  can  be  written  from  the  signal 
relationship. 

Yt  =  X  Hs  (Ht  +  Hb)  +  Nt 

Yb  =  X  Hs  (  0  +  Hb)  +  Nb 

Yc  =  X  Hs  (He  +  Hb)  +  Nc 

where  Yt  =  signal  out  of  network  analyzer  in 

tcurget  measurement 

Yb  =  signal  out  of  network  analyzer  in 
background  measurement 

Yc  =  signal  out  of  network  analyzer  in 
calibration  sphere  measurement 

X  =  RF  excitation 

Hs  =  frequency  response  of  measurement  system 
dependency 

Hb  =  frequency  response  of  background 

An  estimator  of  the  frequency  response  of  the  target,  Hte» 
can  be  derived  by  a  simple  algebraic  manipulation  of 
Equations  (3) -(5).  This  will  be  termed  the  calibration 
process . 


(3) 

(4) 

(5) 


Hte  =  Hcc  (Yt  -  Yb)  /  (Yc  -  Yb)  (6) 

The  estimator  becomes  exact  as  the  measurement  noise 
approaches  zero. 

As  shown  in  Equaton  (6),  the  target  frequency  response 
can  be  calculated  from  the  three  measured  data  and  the 
con5)uted  frequency  response  of  the  calibration  sphere,  Hcc* 
The  measurement  system  dependency,  Hs,  cancels  out  in 
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deriving  Equation  (6).  From  the  definition  in  Equation  (1), 
the  radar  cross  section,  o,  is  given  by 

a  =  lHt|2  (7) 

Raw  data  taken  from  the  network  analyzer,  Yt#  is  shown  in 
Figvire  10.  The  cuitenna  coupling  is  substantially  larger  than 
the  target  scattering  at  the  lower  frequencies  so  that  the 
backgro\ind  signal  dominates  over  the  target  signal.  Figure 
11  shows  the  target  scattering  after  the  background  is 
siibtracted  out;  the  cuitenna  coupling  and  chamber  clutter  are 
eliminated.  The  measurement  system  dependency  is  removed  in 
the  calibration  process  of  Equation  (6),  and  the  result  is 
shown  in  Figure  12.  The  spikes  in  the  data  are  caused  by  a 
residual  antenna  coupling  and  chamber  clutter,  and  will  be 
greatly  reduced  by  time-domain  gating  to  be  discussed  in  the 
following  section. 

C .  DIGITAL  SIGNAL  PROCESSING 
1 .  Window  Shape 

The  frequency-domain  data  are  weighted  with  a  window 
fxinction  prior  to  inverse  FFT  in  order  to  reduce  sidelobes  in 
the  time  domain.  Among  many  window  shapes  available  in  the 
MatLab  signal  processing  library,  the  Hamming  window  [Ref. 
6],  which  is  a  raised  cosine  function,  was  determined  to  be 
best  suited  because  the  effects  of  windowing  can  be  removed 
by  multiplying  the  data  with  an  inverse  of  the  window 
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function.  Other  windows  have  zero  end  values  so  that  the 
inverse  can  not  be  taken.  The  actual  equation  of  the  Hainming 
window  in  the  library  is  given  in  Equation  (8) . 

W  =  0.54  -  0.46*cos{2*pi*n/ (N-l> )  (8) 

where  w  =  weighting  factor  in  the  window  function 
N  =  total  number  of  sample  points 
n  =  index  ranging  from  0  to  N-1 

The  data  in  Figure  12  are  Hamming -windowed  and  shown 
in  Figure  13. 

2 .  Zero-padding  and  Complex  Conjugate 

After  the  data  are  weighted  with  the  Hamming  window 
function,  100  data  points  with  zero  value  are  inserted  in 
front  of  the  windowed  data  to  fill  in  the  missing  dc-2ghz 
data,  and  123  zeros  are  appended  to  the  end  of  the  data  to 
make  the  total  number  of  data  points  1024,  thus  allowing  use 
of  radix-2  FFT  algorithm.  The  complex  conjugates  of  the  1024 
data  points  are  added,  making  a  total  of  2048  data  points. 
The  complex  conjugate  is  necessary  to  make  the  transformed 
time-domain  data  a  time  transient  response.  The  layout  of 
the  2048  data  points  is  illustrated  in  Figure  14. 

3 .  Fourier  Transforas#  Gating  and  Unwindowing 

The  2048  data  points  in  the  frequency  domain  are 
transformed  to  the  time  domain  with  an  IFFT  routine  in  the 
MatLab  library.  The  routine,  IFFT(x),  is  the  inverse 
discrete  Fourier  transform  of  vector  X.  The  IFFT  employs  a 
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Frequency  Response  after  Windowing 
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fast  radix-2  algorithm  when  the  length  of  X  is  an  integer 
power  of  two,  and  a  slower  algorithm  when  the  length  of  X  is 
not  a  power  of  two.  The  IFFT  operation  performed  on  the  2048 
complex  data  points  yields  another  set  of  2048  complex 
numbers  of  which  the  imaginary  part  is  zero,  and  the  real 
part  represents  the  time  transient  response. 

In  the  time-domain  plot,  the  position  where  the 
target  appears  depends  on  the  phase  reference  in  the 
confutation  of  the  calibration  sphere.  In  order  to  place  the 
antenna  coupling  at  time  zero  and  thus  make  the  time 
transient  response  plot  easier  to  read,  the  data  are  shifted 
589  cells  as  shown  in  Figure  15.  The  wrap-around  nature  of 
the  Fourier  transform  output  permits  the  shifting  of  data 
without  any  loss  or  discontinuity.  Shifting  data  in  the  time 
domain  is  equivalent  to  multiplication  in  the  frequency 
domain  by  a  linear  phase  term  and  does  not  affect  the 
magnitudes . 

After  shifting,  the  time  domain  data  are  gated  to 
remove  residuals  of  antenna  coupling  and  any  return  from  the 
chamber  which  has  not  been  eliminated  completely  in  the 
background  subtraction.  The  gate  is  centered  at  the  600th 
cell,  and  the  width  is  400  cells.  The  gate  is  centered  on 
the  target,  and  the  gate  width  corresponds  to  approximately 
10  nanoseconds. 

The  gated  data  are  transformed  back  to  the  frequency 
domain  with  an  FFT  routine  in  the  MatLab  library.  Because 
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Figure  15.  Time  Response  after  Shifting 


zeros  and  conjugates  have  been  appended  to  the  frequency 
data,  the  801  data  points  between  the  101st  cell  and  the 
901st  cell  represent  the  2-18  GHz  frequency  range.  Finally, 
an  inverse  of  the  Hamming  window  is  applied  to  the  801  data 
points  to  remove  the  effect  of  windowing.  The  final  data  are 
shown  in  Figure  16  and  the  sharp  spikes  in  Figure  12  are 
eliminated. 
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Frequency  Response  after  Unwlndowlng 


V .  SYSTEM  VALIDATION 

A .  CANONICAL  TARGETS 

A  set  of  metal  spheres  was  measured  in  the  NPS  TESL 
anechoic  chamber  in  order  to  validate  the  fidelity  of  the 
frequency-domain  scattering  measurement  system.  Metal 
spheres  were  chosen  as  validation  targets  because;  1)  the 
exact  solution  is  readily  available  from  the  Mie  series 
solution,  2)  the  RCS  is  aspect-independent  so  that  errors  in 
target  positioning  are  negligible,  and  3)  the  time  transient 
response  exhibits  a  delayed  creeping  wave  which  can  be  used 
to  validate  the  time-domain  transformation. 

The  computer  program,  MIE. BAS  written  by  M.  Morgan  for 
the  Mie  series  solution  [Ref.  2],  has  been  modified  and  named 
MIE2.BAS.  Since  the  original  MIE. BAS  calculates  the 
scattered  field  observed  at  the  receiving  antenna,  the  free- 
space  spread  factor  defined  in  Equation  (2)  has  been  added  to 
the  new  MIE2.BAS  program.  The  output  of  MIE2.BAS  is  used  in 
the  post-processing  program,  PROC.M,  for  calibrating  the 
measured  data  by  supplying  the  Hcc  iii  Equation  (6). 

B .  EXPERIMENTAL  RESULTS 

A  series  of  scattering  measurements  was  performed  on 
metal  spheres  in  the  frequency  range  of  2  to  18  GHz.  The 
spheres  were  15.24,  12.07,  7.94,  and  2.54  centimeters (cm)  in 


37 


diameter.  Setup  of  the  HP-8510B  was  in  STEP  mode  sweep,  801 
points,  and  averaging  factor  of  128.  Antennas  were  oriented 


for  horizontal  polarization.  The  12.07  cm  sphere  was 
arbitrarily  chosen  as  a  calibrator  and  the  other  spheres  as 
targets.  A  test  matrix  is  svimmarized  in  Table  2. 

TABLE  2.  TEST  MATRIX  OF  SCATTERING  MEASUREMENTS 


Three  measurement  steps  were  taken  for  each  target: 
calibration,  background,  and  target.  With  a  128  averaging 
factor,  60  seconds  were  required  for  a  complete  sweep  of  801 
frequency  points,  15  seconds  to  transfer  data  from  the  HP- 
8510B  to  the  PC,  and  10  seconds  to  store  this  data  to  the 
hard  disk.  When  no  averaging  is  taken,  sweeping  through  801 
frequency  steps  takes  40  seconds,  most  of  it  is  for  phase¬ 
locking  the  receiver.  As  the  averaging  factor  increases,  the 
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sweep  time  also  increases  progressively  to  the  maximum  of  700 
seconds  for  the  maximum  averaging  factor  of  4096.  Trade-offs 
between  the  averaging  factor  and  the  measurement  time  must  be 
considered  in  selecting  the  averaging  factor. 

The  MIE2.BAS  program  was  executed  for  the  12.07  cm 
calibration  sphere.  The  computation  includes  the  small 
bistatic  angle  between  two  antennas  and  provides  the 
magnitude  and  phase  of  Hcc*  These  data  were  stored  under 
filenames  of  "mag.dat"  and  "faz.dat". 

Results  of  the  "PROC.M"  post-processing  program  are  shown 
in  Figures  17  through  22.  Figures  17,  19,  euid  21  show  the 
frequency  responses  of  the  target  spheres,  with  the 
measurement  data  plotted  in  dotted  lines,  and  theoretical 
values  from  the  Mie  solution  in  solid  lines.  Comparison  of 
the  frequency  responses  shows  good  agreement,  and  thus 
confirms  the  validity  of  the  measurement  system.  Figures  18, 
20,  and  21  show  the  corresponding  time  responses.  The  time 
responses  represent  the  scattered  waveform  resulting  from  the 
incident  impulse  formed  by  the  inverse  transform  of  the  2-18 
GHz  bandlimited  signal  weighted  by  the  Hanuning  window.  The 
time  responses  show  the  main  specular  return,  from  the  front 
face  of  the  sphere,  followed  by  the  creeping -wave,  delayed  in 
time  corresponding  to  a  single  circ\imnavigation. 

A  detailed  comparison  between  the  measurement  data  and 
computational  data  in  the  frequency  responses  indicates  a 
minor  discrepancy  of  less  than  1  dB  in  magnitude.  This 
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Figure  1.8.  Time  Transient  Response  of  a  2.54  cm  Diameter  Sphere 
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figure  19.  frequency  Response  of  a  7.94  cm  Diameter  Sphere 
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frequency  Response  of  a  15.24  cm  Diameter  Sphere 


Time  Transient  Response 
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Figure  22.  Time  Transient  Response  of  a  15,24  cm  Diameter  Sphere 


discrepancy  remained  even  when  the  averaging  factor  was 
increased  substantially,  showing  that  the  error  results  from 
coherent  sources.  Possible  error  sources  are  1)  interaction 
between  the  sphere  and  the  chamber,  2)  disturbance  of  fields 
on  the  sphere  caused  by  the  target  support  column  touching 
the  sphere,  3)  non-uniform  field  distribution  in  the  target 
zone  caused  by  the  antenna  beam  pattern,  and  4)  imperfections 
in  the  sphere's  surface  radius,  as  well  as  scratches  and 
dents.  Despite  these  possible  error  sources,  the  measurement 
results  are  in  good  agreement  with  the  computational  data, 
thus  confirming  the  validity  of  the  measurement  system  and 
post-processing  algorithm. 
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VI.  CONCLUSIONS 


A .  SUNNART 

A  new  frequency-domain  electromagnetic  scattering 
measurement  capability  has  been  developed  for  the  NFS  TESL 
cunechoic  chamber.  The  new  system  takes  advantage  of  the  HP- 
8510B  network  analyzer  for  its  broad  frequency  coverage,  wide 
dynamic  rcuige,  and  automated  control  of  sub-systems.  A  pair 
of  AEL  H1498  horns  operating  at  2  to  18  GHz  is  installed  as 
transmitting  and  receiving  antennas.  An  HP-823 OOC  BASIC 
Language  Processor  installed  on  a  COMPAQ  Deskpro-386  PC 
enables  remote  control  of  the  HP-8510B  and  data  acquisition 
over  the  HPIB  bus.  A  post -processing  algorithm  has  been 
developed  in  MatLab  for  backgroxand  subtraction,  calibration, 
windowing,  adding  zeros  and  conjugates,  IFFT,  shifting, 
gating,  and  FFT. 

The  new  system  was  utilized  to  measure  the  scattered 
fields  from  several  different-size  spheres.  Good  agreement 
between  the  measurements  and  computational  Mie  series 
demonstrated  the  validity  of  the  measurement  system  and  the 
post-processing  algorithm. 

B .  FUTURE  CONSIDERATIONS 

The  antenna  component  needs  to  be  improved  to  enhance  the 
measurement  fidelity  and  expand  the  frequency  bandwidth.  For 


low  frequencies,  the  increased  antenna  coupling  limits  the 
dynamic  range  and  SNR  of  the  system.  The  coupling  between 
antennas  must  be  reduced  by  either  moving  them  apart, 
inserting  RAM  between  them,  finding  a  new  antenna  that 
controls  the  coupling  better,  or  gating  the  RF  in  hardware. 
The  antenna  also  limits  the  frequency  bandwidth.  The 
frequency  coverage  of  the  current  system  can  be  increased  to 
26  GHz  by  upgrading  the  antennas  and  the  directional  coupler. 

When  the  new  in^roved  antennas  are  in5>lemented  to  reduce 
the  antenna  coupling,  the  RAMP  mode  of  the  HP8510B  may  be 
used  to  speed  up  the  measurement  time.  In  the  RAMP  mode,  the 
network  analyzer  sweeps  through  the  801  frequency  points  in  a 
fraction  of  a  second. 

One  of  the  advantages  of  the  HP-8510B  based  system  is  its 
flexibility  for  expansion  into  the  mil lime ter -wave  region. 
With  a  simple  addition  of  a  second  RF  source,  frequency 
multipliers,  and  millimeter-wave  components  such  as  a 
directional  coupler  and  antennas,  the  frequency  range  of  the 
RCS  measurement  can  be  extended  to  110  GHz. 

One  last  consideration  for  the  future  is  incorporating  a 
pedestal  for  the  rotation  of  targets.  The  ability  to  rotate 
the  target  in  precise  angle  increments  and  correcting 
coherent  data  at  different  aspect  angles  allows  the 
generation  of  two-dimensional  scattering  images  of  the  target 
using  inverse  synthetic  aperture  radar  (ISAR)  techniques. 
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APPENDIX  A 


SOUECB  CODE  OP  THE  “GET_DATA"  PROGRAM 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

ISO 

160 

170 

180 

190 

200 

210 

220 


i  GET_DATA:  THIS  HP  BASIC  PROGRAM  CONTROLS  HP8510 

!  REMOTELY  AND  WRITES  MAG  &  PHASE  DATA  ON  PC 

!  HARO  DISK  UNDER  C:\BLPSUBDIRECTORATE. 

1  WRITTEN  BY  KEN  OH 

!  REVISED  ON  9/2/93 

OPTION  BASE  1 
COM  REAL  Exp_tbl (0:255) 

MASS  STORAGE  IS  "XBLPtDOS^C" 

PRINTER  IS  1 
PRINT  CHR$(12) 

ASSIGN  eHplb  TO  7 
ASSIGN  §Nwa  TO  716 
REMOTE  eHplb 
CLEAR  §MWa 
Exp_tbl(0)=2*(-15) 

FOR  1*0  TO  126 

Exp_tbl ( I+l ) *Exp_tbl ( I ) +Exp_tbl ( I ) 

NEiXT  X 

Exp  tbl(128)*2*(-143) 

P0R”I»128  to  254 

Exp  tbl(I+l)*Exp  tbl(I)+Exp_tbl(I) 

NEXT  I 


230  Main  menu:  i 
240  PRINT  CHR$(12) 

250  OUTPUT  eNwa;"STAR;OUTPACTI;" 

260  ENTER  fNva ; Steu:_f req 

270  OUTPUT  §Nwa;"STOP;OUTPACTI;" 

280  ENTER  iNva ; Stop_f req 

290  OUTPUT  eNwa ; "POIK ; 7UTPACTI ; " 

300  ENTER  eNwa;Npoint 

310  OUTPUT  §Nwa;"POWE;OUTPACTI;" 

320  ENTER  iNWa; Power 
330  OUTPUT  eNwa ; "SWEM? ; " 

340  ENTER  eNwa;Swe$ 

350  OUTPUT  eNwa ; "CALS? ; " 

360  ENTER  eNWa;Calset 

370  OUTPUT  eNWa ;  "AVERON;OUTPACTI ; " 

380  ENTER  eNwa;Avn 
390  OUTPUT  eNwa;"ENTO;" 

400  DISP 

410  PRINT  "Main  Menu" 

420  PRINT  "  " 

430  PRINT  "START  FREQ  IS:";Star_freq/l.E+9;"GH2" 

440  PRINT  "STOP  FREQ  IS:";Stop_freq/l.E+9;"GHz" 

450  PRINT  "NUMBER  OF  POINTS  IS:";Npoint 
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460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 


PRINT  "SOORCB  PONESt  IS: "jINTCPower-i-. 001)  ; 

PRINT  "SHEEP  MODE  IS:";Sve$ 

PRINT  "ACTIVE  CAL  SET  IS:";Cal8et 
PRINT  "AVERAGING  FACTOR  IS:";Avn 
OFF  KEY 

LABEL  "Continu. Sweep"  GOTO  Cont_sweep 
LABEL  "Single  Sweep"  GOTO  Sing_sweep 
LABEL  "Hold  Sweep"  GOTO  Hold_sweep 
LABEL  "Change  #Average"  GOTO  Average 
LABEL  "Record  Data"  GOTO  Call_record 
LABEL  "PlaybackData"  GOTO  Call_playback 
LABEL  "Change  Setup"  GOTO  Setup 
ON  KEY  8  LABEL  "Default  Setup"  GOTO  Default 
GOTO  510 
! 
i 


ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 
ON  KEY 


1 

2 

3 

4 

5 

6 
7 


620  Cont_8weep:  ! 

630  DISP  "CONTINUOUS  SWEEP" 

640  OUTPUT  SNwa ; "CONT ; " 

650  GOTO  Main__iiienu 
660  Sing_8weep:  ! 

670  DISP  "SINGLE  SWEEP  _  PLEASE  WAIT" 

680  OUTPUT  eNWa ; "SING ; " 

690  GOTO  Main_nenu 
700  Hold_8weep:  1 
710  DISP  "HOLD  SWEEP" 

720  OUTPUT  §Nwa;"HOLD;" 

730  GOTO  Nainjnenu 
740  Average:  ! 

750  INPUT  "ENTER  THE  NUMBER  OF  AVERAGE: " ,Avn 
760  OUTPUT  §Nwa;"AVERON";Avn 
770  GOTO  Nainjnenu 
780  Call_record:  ! 

790  DISP  "RECORDING  DATA -  PLEASE  WAIT." 

800  CALL  Record_data(Npoint} 

810  GOTO  Nainjnenu 

820  Call_playback:  ! 

830  DISP  "READING  DATA _ PLEASE  WAIT." 

840  CALL  Playback__data(Npoint) 

850  GOTO  Main_menu 

860  Setup:  ! 

870  PRINT  CHR$(12) 

880  OUTPUT  §Nwa;"STAR;OOTPACTI;" 

890  ENTER  9Nwa ; Star_f req 

900  OUTPUT  iNwa;"STOP;ODTPACTI;" 

910  ENTER  fNwa;Stop_freq 

920  OUTPUT  §Nwa;"POIN;OUTPACTI;" 

930  ENTER  eNwa;Npoint 

940  OUTPUT  eNwa;"POWE;OUTPACTI;" 

950  ENTER  8Nwa; Power 
960  OUTPUT  eNwa;"SWEM?;" 


50 


970  ENTER  9Nwa;Sv«$ 

980  OUTPUT  §NWa;"CALS?;" 

990  ENTER  8NVa;Cal8et 
1000  PRINT  "Setup  Menu" 

1010  PRINT  "  " 

1020  PRINT  "START  FREQ  IS:";Star_£req/l.E+9;"GH2" 

1030  PRINT  "STOP  FREQ  IS:";Stop_freq/l.E+9;"GHz" 

1040  PRINT  "NUMBER  OF  POINTS  IS:";Npoint 

1050  PRINT  "SOURCE  POWER  IS:";INT(Power+.001)  ;"dB]n" 

1060  PRINT  "SWEEP  MODE  IS:";Sve$ 

1070  PRINT  "ACTIVE  CAL  SET  IS:";Calset 

1080  OFF  KEY 

1090  ON  KEY  1  LABEL  "Save  State"  GOTO  Save_state 

1100  ON  KEY  2  LABEL  "Recall  State"  GOTO  Recall_state 

1110  ON  KEY  3  LABEL  "Change  Freg"  GOTO  Freg_range 
1120  ON  KEY  4  LABEL  "No  of  Point"  GOTO  Nojpoint 

1130  ON  KEY  5  LABEL  "Power  Level"  GOTO  Power 

1140  ON  KEY  6  LABEL  "Sweep  Mode"  GOTO  Sweepjnode 
1150  ON  KEY  7  LABEL  "Calibra-tion"  GOTO  Call_cal 
1160  ON  KEY  8  LABEL  "Go  to  MainMenu"  GOTO  Main  menu 
1170  GOTO  1090  ■" 

1180  Recall_8tate: ! 

1190  INPUT  "ENTER  RECALL  STATE  NUMBER: ", Sset 
1200  OUTPUT  §Nwa;"RECA"&VAL$(Sset) 

1210  GOTO  Setup 
1220  Save__8tate : ! 

1230  INPUT  "ENTER  SAVE  STATE  NUMBER: ",Sset 
1240  OUTPUT  §Nwa;"SAVE"&VAL$(Sset) 

1250  GOTO  Setup 
1260  FreG[_range: ! 

1270  INPUT  "ENTER  START  FREQ  IN  GHz : " , Star_f req 
1280  INPUT  "ENTER  STOP  FREQ  IN  GHz : " , Stop_f req 
1290  OUTPUT  9Nwa ;" STAR" ;Star_freq;" GHz?" 

1300  OUTPUT  §Nwa ; "STOP" ?  Stop_f req ; "GHz ; " 

1310  GOTO  Setup 
1320  No_point: 1 

1330  INPUT  "ENTER  NUMBER  OF  POINT  (51,101,201,401  or  801):" 
,Npoint 

1340  OUTPUT  §Nwa ; "POIN" , Npoint 
1350  GOTO  Setup 
1360  Power: 1 

1370  INPUT  "ENTER  SOURCE  POWER  LEVEL  IN  dBm:", Power 
1380  OUTPUT  8NWa  ? "POWE" , Power 
1390  GOTO  Setup 
1400  Sweepjmode:! 

1410  INPUT  "ENTER  1  FOR  STEP  MODE  AND  2  FOR  RAMP  MODE:",Swe 
ep 

1420  IF  Sweep«2  THEN 
1430  GOTO  RampjBode 
1440  ELSE 

1450  GOTO  Step_mode 
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1460  BMD  Zr 
1470  8t^  BOd*:! 

1480  ODlfDT  tllva;"STEP;" 

1490  GOTO  8«tup 
1500  RaapjBod*:! 

1510  OUTPUT  CNva ;  "RAMP;  •* 

1520  GOTO  Setup 
1530  Call_cal:! 

1540  INPUT  "ENTER  CAL  SET  NO  (1>8) ,  OR  0  FOR  NEW  CAL,  OR  9 

FOR  CAL  OFF",Cn 

1550  IF  Cn»0  THEN 

1560  CALL  Calibrate (Npoint) 

1570  GOTO  Maim 

1580  END  IF 

1590  IF  Cn<9  THEN 

1600  OUTPUT  fNwa;"CORRON;" 

1610  OUTPUT  6Nwa;"CALS"4VAL$(Cn) 

1620  ELSE 

1630  OUTPUT  §Nwa ; "CORROFF ; " 

1640  END  IF 
1650  GOTO  Setup 
1660  Default: ! 

1670  DISP  "PRESETTING  THE  NWA  ...  PLEASE  WAIT." 

1680  OUTPUT  §Nwa;"PRES;" 

1690  OUTPtrr  iNwa;"STAR  2  GHz;STOP  18  GHz;" 

1700  OUTPUT  eNwa;"POIN801;" 

1710  OUTPUT  §Nwa;"STEP;" 

1720  OUTPUT  «NWa;"CONT;" 

1730  OUTPUT  «Nwa;"OUTPERRO;" 

1740  GOTO  Main  nenu 
1750  1 ! 

1760  STOP 
1770  END 
1780  !! 

1790  SUB  Calibrate (Npoint) 

1800  OPTION  BASE  1 

1810  ASSIGN  SNwa  TO  716 

1820  ASSIGN  8Nwa_data2  TO  716;FORMAT  OFF 

1830  INTEGER  Preamble, Size, Cals 

1840  DIM  Rooa(801,2) ,Reference(801,2) ,Source(801,2) 

1850  DISP  CHR$(129)&"MEASURE  ROOM, "&CHR$ (128) &"  THEN  PRES 

S  CONTINUE" 

I860  BEEP 

1870  LOCAL  716 

X880  PAUSE 

1890  OUTPUT  6Nwa;"FORM3;OUTPRAHl" 

1900  ENTER  6Nwa_data2 ; Preeuoble ,  Size ,  Room  ( * ) 

1910  1 

1920  DISP  CHR$(129)&"MEASURE  REFERENCE, "&CHR$ (128) &"  THEN 

PRESS  CONTINUE" 

1930  BEEP 
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1940  LOCAL  716 

1950  PAUSE 

I960  OUTPUT  tNwa;"FORM3;OUTPRAWl" 

1970  ENTER  9Nwa_data2 ; Preaable , Size , Reference ( * ) 

1980  1 

1990  DISP  CHR$(131)&"SUBTRACTIN6  OUT''&CHR$(128) &”  THE  EFF 

ECTS  OF  THE  ROOM" 

2000  FOR  M«1  TO  2 

2010  FOR  N«1  TO  Npoint 

2020  Reference  (N,M) -Reference  (M,M)  -*Rooin(N,M) 

2030  NEXT  N 

2040  NEXT  M 

2050  ! 

2060  INPUT  "ENTER  CAL  SET  NUMBER  (1-8)  AND  PRESS  CONTINUE 
" , Cals 

2070  BEEP 

2080  IF  Cals«0  THEN  GOTO  2060 

2090  DISP  CHR$(131)&"LOADIN6"&CHR$(128)&"  CAL  COOEFFICIEN 

TS" 

2100  OUTPUT  9Nwa;"CALl;CALISlll" 

2110  OUTPUT  iNwa;"FORM3;INPUCALC01" 

2120  OUTPUT  §Nwa_data2 ; Preamble, Size, Room(*) 

2130  OUTPUT  9Nva;"FORM3;INPUCALC02" 

2140  OUTPUT  @Nwa_data2;Preamble, Size, Source (*) 

2 150  OUTPOT  8Nwa ; "FORM3 ; INPUCALC03 " 

2 160  OUTPUT  9Nwa_data2 ; Preamble , Size , Reference ( * ) 

2170  OUTPUT  eNwa; " SAVC; CALS" &VAL$ (Cals) 

2180  DISP 

2190  LOCAL  716 

2200  SUBEND 
2210  !! 

2220  SUB  Record_data (Npoint) 

2230  OPTION  BASE  1 

2240  ASSIGN  iHpib  TO  7 

2250  ASSIGN  9Nwa  TO  716 

2260  PRINTER  IS  1 

2270  REAL  Freq(l:801) ,Mag(l:801) ,Phase(l:801) 

2280  REAL  Block(l:801, 1:2) 

2290  INTEGER  Binary(l:801,0:2) 

2300  INTEGER  I, Points, M,N 

2310  COM  Exp__tbl(*) 

2320  REMOTE  8Hpib 

2330  ABORT  SHpib 

2340  CLEAR  eNwa 

2350  Points^^Npoint 

2360  Enter_data(8Nwa, Points, Binary (*) ) 

2370  Convertb( Points, Binary (*) ,Block(*) ) 

2380  Transform (Points, Block (*) ,Mag(*) ,Phase(*) ) 

2390  CLEAR  8Nwa 

2400  BEEP 

2410  PRINT  CHR$(12) 
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2430  DIBB 

2430  GINIT  I *******#*PL0T  GRAPH 

2440  GCLSAR 

2450  GRAPHICS  ON 

2460  VIEHPORT  10,140,60,100 

3470  FRAME 

2480  Maxv— 10000 

2490  N-Npolnt 

2500  FOR  1-1  TO  N 

2510  Naxv>NAX(Maxv,Mag(I)) 

2520  NEXT  I 

2530  WINDOW  l,N,Maxv-40,Maxv-t'5 
2540  FOR  I-l  TO  N 

2550  MOVE  1,-200 

2560  DRAW  I,Mag(I) 

2570  NEXT  I 

2580  VIEWPORT  10,140,10,50 
2590  FRAME 

2600  WINDOW  1,N,-180,180 

2610  FOR  I-l  TO  N 

2620  MOVE  1,0 

2630  DRAW  I, Phase (I) 

2640  NEXT  I 

2650  1 

2660  PAUSE 

2670  GCLEAR 

2680  1 

2690  IFILE  CREATOR 

2700  INPUT  "ENTER  NEW  FILE  NUMBER  =",Istrt 

2710  Filen$-"" 

2720  DISP  "WRITING  DATA  TO  DISK. . . .  PLEASE  WAIT" 

2730  Filen$«"F"&VAL$(Istrt)S".DAT" 

2740  CREATE  BDAT  Filen$, 26, 256 

2750  ASSIGN  fPath  TO  Filen$ 

2760  OUTPUT  §Path;Mag(*) ,Phase(*) 

2770  ASSIGN  §Path  TO  * 

2780  Filen$="D"&VAL$(Istrt)6".DAT" 

2790  CREATE  Filen$,256 

2800  ASSIGN  ePath2  TO  Filen$ 

2810  FOR  I-l  TO  N 

2820  OUTPUT  §Path2  USING  "K";I,CHR$(44) ,Mag(I) ,CHR$(44) 

, Phase  (I) 

2830  NEXT  I 

2840  GCLEAR 

2850  ASSIGN  §Path2  TO  * 

2860  SUBEND 
2870  !! 

2880  DEF  FNAtan(X,Y,P' 

2890  IF  X-0  THEN  X»>^.E-9 

2900  D-ATN(Y/X) 

2910  IF  (X<0)  AND  (Y>0)  THEN  P-D+PI 
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2920  IP  (X<0)  AND  (Y<0)  THEN  P-D-PI 

2930  IF  X>0  THEN  P»D 

2940  RETURN  P 

2950  FNENO 

2960  !! 

2970  SUB  Print_data( INTEGER  N,REAL  X(*),REAL  Y(*),REAL  H(*) 

) 

2980  INTEGER  I 

2990  FOR  I»1  TO  128 

3000  PRINT  USING  "4 (5D.DD) ";I,X(I) , Y(I) ,M(I) 

3010  NEXT  I 

3020  SUBEND 
3030  !! 

3040  SUB  Enter_data(eN, INTEGER  N,D(*)) 

3050  INTEGER  Preanble , Bytes 

3060  OUTPUT  fN;"FORMl;  OUTPDATA;" 

3070  ASSIGN  §N;FORMAT  OFF 

3080  ENTER  §N;Preeunble; Bytes 

3090  N>Bytes/6 

3100  REDIM  D(l:N,0:2) 

3110  ENTER  eN;D(*) 

3120  ASSIGN  iN; FORMAT  ON 

3130  SUBEND 
3140  !  ! 

3150  SUB  Transform (INTEGER  N,REAL  D(*) ,M(*) ,P(*) ) 

3160  INTEGER  I 

3170  REAL  R,Im,Mag 

3180  DEG 

3190  FOR  1=1  TO  N 

3200  R=D(I,1) 

3210  Im=D(I,2) 

3220  Mag=SQR(R*R+Im*Im) 

3230  M(I)=20*LGT(Mag) 

3240  P(I)=2*ATN(Im/(R+Mag)) 

3250  NEXT  I 

3260  SUBEND 
3270  !! 

3280  SUB  Convertb( INTEGER  N,B(*),REAL  D(*)) 

3290  INTEGER  I 

3300  REAL  Exp 

3310  COM  Exp_tbl(*) 

3320  FOR  1=1  TO  N 

3330  Exp=Exp_tbl(BINAND(B(I,2) ,255)} 

3340  D(I,l)=B(I,l)*Exp 

3350  D(I,2)=B(I,0)*Exp 

3360  NEXT  I 

3370  SUBEND 
3380  !! 

3390  DEF  FNLgt(X,Y) 

3400  Z=X“2+Y*2 

3410  IF  Z>. 000001  THEN 
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3420  VlO9»10*L0T(Z) 

3430  ELSE 

3440  Vlog— 200 

3450  END  IP 

3460  EETUKN  Vlog 

3470  FNEND 
3480  11 

3490  SDB  Playbac)c__data  (Npoint) 

3500  OPTION  BASE  1 

3510  ASSIGN  eHpib  TO  7 

3520  ASSIGN  ONWa  TO  716 

3530  PRINTER  IS  1 

3540  REAL  Freq(l:801) ,Mag(l:801) ,Phase(l:801) 
3550  REAL  Bloc]c(l:801,l:2) 

3560  INTEGER  Binary (1:801, 0:2} 

3570  INTEGER  I, Points, M,N 

3580  cam  Exp_tbl(*) 

3590  INPUT  "ENTER  FILE  NUMBER 

3600  FnaBe$»"F"&VAL$(M\m) 

3610  OUTPUT  KBD;Claar$; 

3620  OUTPUT  KBD;Hoiiia$; 

3630  ASSIGN  8Path  TO  Pname$ 

3640  ENTER  SPath ; Nag (*), Phase (*} 

3650  PRINT  CHR$(12) 

3660  DISP 

3670  GINIT  1 *********PLOT  GRAPH 

3680  GCLEAR 

3690  GRAPHICS  ON 

3700  VIEWPORT  10,140,60,100 

3710  FRAME 

3720  Maxv»-10000 

3730  N-Npolnt 

3740  FOR  1*1  TO  N 

3750  Naxva4(AX(Maxv,Nag(I) ) 

3760  NEXT  I 

3770  WINDOW  l,N,Maxv-40,Maxv+5 

3780  FOR  1*1  TO  N 

3790  MOVE  1,-200 

3800  DRAW  I,Mag(I) 

3810  NEXT  I 

3820  VIEWPORT  10,140,10,50 

3830  FRAME 

3840  WINDOW  1,N,-180,180 

3850  FOR  1*1  TO  N 

3860  MOVE  1,0 

3870  DRAW  I, Phase (I) 

3880  NEXT  I 

3890  PAUSE 

3900  GCLEAR 

3910  SUBEND 
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APPENDIX  B . 


SOURCE  CODE  OF  THE  "PROC.M"  PROGRAM 


%  Post-processing  program 
%  REV.  6/24/93  Ken  Oh 

clear,  clg 

filenaBe«input( 'Enter  Cal  File:  dlOl 

eval(['load  ' ,£ilename, ' .dat* ] ) 
datal»eval (filename) ; 

£ilename«input( 'Enter  BG  File:  dl02  ','s'); 

eval(['load  ', filename, ' .dat* ] ) 

data2«eval(filen^uBe) ; 

filenamesinput( 'Enter  Target  File:  dl03  ','s'); 
eval ( [ ' load  ' , filename, ' . dat ' ] ) 
data3»eval( filename) ; 

filenameaBinput( 'Enter  Computed  Sphere  Mag  File:  mag  ','s') 
eval(['load  ', filename, ' .dat' ] ) 
data4»eval (filename) ; 

filename«input( 'Enter  Computed  Sphere  Phase  File:faz  ','s'} 
eval(C'load  ' ,£ilename, ' .dat' ] ) 
dataS-eval (filename) ; 

filename«input( 'Enter  Target  Mie  File:  ml,  m2,  or  m3  ','s') 
eval(['load  ' ,filename, ' .dat' ] ) 
datae^eval ( f i lename ) ; 


j“sqrt(-l);,  npoint=801; 
freq»[2: (18-2) /800:18] ' ; 
t-C0:50/800:50] '; 
t2048»[0:50/2047:50] '  ; 
tzoom-[ 10: 10/409:20] ' ; 


%  ******************************* 
ml»zeros(l,npoint) ; 
plszeros ( 1 , npoint ) ; 
for  i>l: npoint 

ml ( i) «datal ( i , 2 ) ; 
pl(i)*datal(i,3) ; 

end; 

ml«10. . “ (ml/20) ; 
pl«pl*pi/180; 

dl=ml.*cos(pl)+j*ml.*sin(pl) ; 


Calibration  Data 


%  in  db 
%  in  degree 

%  in  linear  unit 
%  in  radian 
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Background  Data 


%  ********4l^***************** 

■2«z«ros ( 1 , npolnt) ; 
p2«zaros(l,i^lnt) ; 
for  i*l: npolnt 

»2(i)»data2(i,2) ; 
p2(i)«data2(l,3) ; 

and; 

a2»10. .^(B2/20) ; 
p2*p2*pi/180; 

d2-<B2.*cos(p2)+j*n2.*sin(p2) ; 


%  in  db 
%  in  degree 

%  in  linear  unit 
%  in  radian 


%  ******************************* 

■3azeros ( 1 , npoint ) ; 
p3azero8(l,npoint) ; 
for  ial: npoint 

B3(i)«data3(i,2) ; 
p3  ^)«>data3  (i,3) ; 

end; 

B3-10. .~(m3/20) ; 
p3«p3*pi/i80; 

d3»®3.*cos(p3)+j*m3.*sin(p3) ; 


Target  Data 


%  in  db 
%  in  degree 

%  in  linear  unit 
%  in  radian 


%  *******************************  Computed  Sphere  Data 
B4»zaros ( 1 , npoint) ; 
p4»zero8 ( 1 , npoint) ; 
for  i»l: npoint 

b4  (i) »data4  ( i-i-4 ) ; 
p4  ( i)  »>data5  ( 1+4) ; 

end; 

d4«qii4 .  *cos  (p4 )  +  j *m4 .  *sin (p4 ) ; 

%  *******************************  Target's  Mie  Solution 
B6«zeros(l, npoint) ; 
for  i»l: npoint 

m6(i)«data6(i+4) ; 

end; 


clg 

plot(freq,B3) 

title ('Target  Raw  Data') 

xlabel ( ' frequency  (GHz) ' ) ,  ylabel ( 'magnitude  (volts) ' ) 

tprintlj 

pause 

%  ******************************* 

dl2»dl**d2;  %  Subtracting  BG  from  Cal 

d32Bd3~d2;  %  Subtracting  BG  from  Target 

clg 

plot(freq,abs(d32) ) 
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title (' Frequency  Response  after  Background  Subtraction') 
xlabel (' frequency  (GHz)'),  ylabel( 'magnitude  (volts)') 
%printlj 
pause 

%  *******************************  Calibration 
f32c»d32./dl2.*d4; 

clg 

plot(freq,ab8(f32c) ) 

title (' Frequency  Response  after  Calibration') 

xlabel ( ' frequency  (GHz) ' ) ,  ylabel ( 'magnitude  (volts) ' ) 

%printlj 

pause 

%  *******************************  Hamming  window 
win»haffiming(npoint) ; 
f32cw»f32c.*win' ; 

clg 

plot  ( f  r  eq ,  win) 

title  ( '  Hzunming  Window ' ) 

xlabel ( ' frequency  (GHz) ' ) ,  ylabel ( 'weighting  factor ' ) 

%printlj 

pause 

clg 

plot(freq,abs(f32cw) ) 

title ( ' Frequency  Response  after  Windowing ' ) 

xlabel ( ' frequency  (GHz) ' ) ,  ylabel ( 'magnitude  (volts) ' ) 

tprintlj 

pause 


%  ****************************  Adding  zeros  and  conjugates 
ff^zeros (1:2048) ; 
f f ( [ 101 : 901] ) =f 32cw; 

ff( [1149: 1949] )=conj (f32cw(801:-l:l) ) ; 
clg 

plot(abs(ff ) ) 

title (' Frequency  Response  after  Adding  zeros  and  Conjugates' 

) 

xlabel (' index  number ') ,  ylabel ( 'magnitude  (volts)') 

%printlj 

pause 

tt=ifft(ff);  %  Inverse  FFT  to  time  domain 

clg 

plot (t2048, real (tt) ) 

title ('Time  Response  before  Shifting') 

xlabel ('time  (nsec)'),  ylabel ( 'magnitude  (volts)') 
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tprintlj 

pauM 


t«iqp->z«ros  (1:1459) ; 
t«ip>tt(  [1:1459]) ; 
tt([l:589])-tt([1460:2048]) ; 

%  Shifting  data  to  put  cross-talk  at  tine  zero, 
tt ([ 590 : 2048 ]) *tenp ; 
clg 

plot (t2048, real (tt) ) 

title ('Tine  Response  after  Shifting') 

xlabel ( 'tine  (nsec) ' ) ,  ylabel ( 'nagnitude  (volts) ' ) 

%printlj 

pause 

gl«zero8 (1:400) ; 
g2»zeros (1:1248) ; 

tt( [1:400] )»gl;  %  Gate  width  «  400  cells 

tt( [801:2048] )»g2; 

clg 

plot(t2048,real(tt) ) 

title ('Tine  Response  after  Gating') 

xlabel ('tine  (nsec)'),  ylabel ( 'magnitude  (volts)') 

tprintlj 

pause 

clg 

plot (tzoom, real (tt( [411:820] ) ) ) ; 

title( 'Zoom-in  View  of  Time  Response') 

xlabel ( 'time  (nsec) ' ) ,  ylabel ( 'magnitude  (volts) ' ) 

%printlj 

pause 

ffg*fft(tt);  %  FFT  to  frequency  response 

clg 

plot(abs(ffg) ) 

title (' Frequency  Response  Showing  All  2048  points') 
xlabel ( ' index  number ' ) ,  ylabel ( 'magnitude  (volts) ' ) 

%printlj 

pause 

f32cg»zeros (1:801) ; 

f32cg«‘ffg(  [101:901] ) ;  %  throw-away  zeros  and  conjugates 

clg 

plot(freq,abs(f32cg) ) 

title ( ' Frequency  Response  from  2  GHz  to  18  GHz ' ) 
xlabel ( ' frequency  (GHz) ' ) ,  ylabel ( 'magnitude  (volts) ' ) 
tprintlj 
pause 

win__inv*i . .  /win ' ; 
f  32cgsf  32cg.  *win__inv; 
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clg 

plot ( f r«q , win^inv) 

title (' Inverse  of  Hanning  window') 

xlabel ( ' frequency  (GHz) ' ) ,  ylabel ( 'weighting  factor ' ) 

tprintlj 

pause 

clg 

plot(freq,abs(f32cg) ) 

title ( ' Frequency  Response  after  Unwindowing ' ) 

xlabel ( ' frequency  (GHz) ' ) ,  ylabel ( 'nagnitude  (volts) ' ) 

%printlj 

pause 

clg 

axis([2  18  -50  -10]) 

plot(freq,20*logl0(abs(f32cg) ) , ' — ' ,freq,20*logi0(n6) , '- 
title ( ' Frequency  Response ' ) 

xlabel ( ' frequency  (GHz) ' ) ,  ylabel ( 'RCS  (dBsn) ' ) 

grid 

tprintlj 

axis; 

end; 


61 


LIST  OF  RBFBRBNCBS 


1.  "Operating  and  Progranuaing  Manual:  HP-8510B  Network 
Analyzer",  Hewlett  Packard,  April  1988. 

2.  M.A.  Mariatequi,  "Development,  Calibration,  and 
Evaluation  of  a  Free  Field  Scattering  Range",  Master's 
Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  1983. 

3.  Brent  W.  McDaniel,  "Calibration  and  Evaluation  of  a  Free 
Field  Scattering  Range  Using  Wideband  Pulse 
Anplification" ,  Master's  Thesis,  Naval  Postgraduate 
School,  Monterey,  CA,  1985. 

4.  "Antenna,  Antenna  Masts,  Mounts  and  Adapters",  American 
Electronic  Laboratories  (AEL)  Inc.,  1989. 

5.  "Installing  and  Using  HP-BASIC  in  the  MS-DOS 
Environment",  Order  Number  82301-90013,  Hewlett  Packard, 
June  1989 . 

6.  "386-MATLAB,  Ver.  3.5h",  MathWorks  Inc.,  Jan.  1991. 


62 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library,  Code  52 

Naval  Postgraduate  School 
Monterey,  CA  93943-5101 

3 .  Department  Chairman,  Code  EC 

Department  of  Electrical  and  Con5>uter  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943-5121 

4.  Professor  Michael  A.  Morgan,  Code  EC/Mw 
Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 

Monterey,  CA  93943-5100 

5.  Professor  David  C.  Jenn,  Code  EC/Jn 
Department  of  Electrical  and  Con^puter  Engineering 
Naval  Postgraduate  School 

Monterey,  CA  93943-5100 

6.  Professor  Rama  Janaswamy,  Code  EC/Js 
Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 

Monterey,  CA  93943-5100 

7.  Dr.  Dean  Mensa,  Code  P03923 

Naval  Air  Warfare  Center-Weapons  Division 
Point  Mugu,  CA  93042-5000 

8.  Dr.  Peter  Yu,  Code  P03923 

Naval  Air  Warfare  Center-Weapons  Division 
Point  Mugu,  CA  93042-5000 

9.  Technical  Library,  Code  P6433A 
BLDG  511A 

Naval  Air  Warfare  Center -Weapons  Division 
Point  Mugu,  CA  93042-5000 


2 


2 


1 


4 


1 


1 


2 


1 


1 


63 


